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qq ' Abstract. We obtained N-band observations of the Apollo asteroid 2002 NY40 during its close Earth fly-by in 

August 2002 with TIMMI2 at the ESO 3.6 m telescope. The photometric measurement allowed us to derive a 

£SJ ■ radiometric diameter of 0.28±0.03 km and an albedo of 0.34±0.06 through the near-Earth asteroid thermal model 

\Q ' (NEATM) and a thermophysical model (TPM). The values are in agreement with results from radar data, visual 

and near-IR observations. In this first comparison between these two model approaches we found that the empirical 

NEATM beaming parameter r\ = 1.0 corresponds to a thermal inertia values of about 100 Jm~ 2 s~ 0,5 K _1 for a 

typical range of surface roughness, assuming an equator-on viewing angle. Our TPM analysis indicated that the 

surface of 2002 NY40 consists of rocky material with a thin or no dust regolith. The asteroid very likely has a 

O I 1 prograde sense of rotation with a cold terminator at the time of our observations. Although both model approaches 

can fit the thermal spectra taken at phase angles of 22° and 59°, we did not find a consistent model solution 

that describes all pieces of photometric and spectroscopic data. In addition to the 2002 NY40 analysis, we discuss 

the possibilities to distinguish between different models with only very few photometric and/or spectroscopic 

jrt , measurements spread over a range of phase angles. 
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F~ , 1. Introduction nitude is H — 19.23 ± 0.2 mag, with variations between 

„„„„ -_.._ .. , , 18.9 mag for the lightcurve maximum and 19.8 mag at the 

I he Apollo asteroid 2002 JNY40 was discovered by , / r> n , i irroni c j c 

-r-r^TT^A^i „ mT1 „ „„„„^,„9\ T , , t-, -, , secondary minimum (sec Bowcll ct al. 1989 for a deiini- 

LINEAR 1 MPEC 2002-O17 2 . It passed the Earth at ,. ,,/ „„ , v N ^^ 

_ _„„_ . TT . „ „ . : , ,. , ,, tion of the H-G system). 

0.0035AU, about 1.3 times the distance to the Moon, 

on August 18th, 2002. This event was closely followed by Howdl et aL ( 200 3» Presented radar images which show 

many amateur and professional observers at a large range that 2002 NY40 looks like tw0 spheroidal units joined to- 

of wavelengths from visible to radar. S ether ' but hi S h resolution optical imaging showed no ob- 

Pravec et al. 12001 found from 21 different observing V10us evidence for a binary structure (H. Mathis 3 ). 

sessions a main lightcurve period of Pi = 19.98 ± 0.01 h 0.3 to 4.0 ^m spectroscopy revealed absorption bands 

with indications for a non-principal axis rotation with a at 1 and 2 /xm due to olivine and pyroxene (Rivkin ct 

non-unique but likely second period of P 2 = 18.43±0.01h. al. I2003JI . The striking similarity with LL6 chondrites, a 

The lightcurve shows sharp bends, indicating the presence subset of ordinary chondrites, indicates that 2002 NY40 

of large non-convex shape features. Its mean absolute mag- could be a progenitor of LL6 meteorites. Based on the 

visible spectrum, Rivkin et al. IJ2T)Tl3"| classified 2002 NY40 

Send offprint requests to: T. G. Miiller as Qdasg asteroidj which is connec ted to a moderately 

Based on observations collected at the European Southern , . , ,, , r r, n or . (r/ fri , , , i njwwM , 

Observatory, Chile; ESO, No. 69.C-0152 hlgh &lbed ° ° f 2 ° " 3 ° % ( Ckrk et aL ^ " 

1 Lincoln Near Earth Asteroid Research project 
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First size estimates were in the order of 700 m to 800 m 
diameter 3,4 , based on magnitude estimates from early as- 
trometric measurements. Later on, adaptive optics sys- 
tems lowered this number to an upper limit of 400 m at 
the time of the observation 5 . Radar data lead to a size 
estimate of about 250 x420 m (E. Howell, personal com- 
munication) consistent with radiometric results (Rivkin 
et al. I2003J1 from 2.5 /im observations. They also gave a 
possible albedo range of 0.15-0.25 with a slight preference 
for the higher value. 

The radiometric method of diameter and albedo deter- 
mination from thermal radiation measurements goes back 
to the early 1970s (references are given in e.g., Morrison 
& Lebofskv ll979f) . A widespread version of this technique, 
the Standard Thermal Model (STM), was published by 
Lebofsky et al. I|1986|l and discussed by e.g., Lebofsky & 
Spencer (^ 



The STM uses an empirical beaming pa- 
rameter rj, which was introduced to adjust the sub-solar 
surface temperature and to account for non-isotropic heat 
radiation (Lebofsky et al. H989|l . The STM is strongly con- 
nected to the IRAS minor planet catalogue (Tedesco et al. 
2002J, with more than 2000 diameter and albedo values 
derived from thermal observations. Morrison l]1976p and 
Cruikshank & Jones Ijl977[l applied this technique for the 
first time to near-Earth asteroids (NEA). However, very 
early on Lebofsky et al. I|1979[l and Veeder et al. I|1989f) 
encountered problems when applying the STM to NEAs. 
They found in several cases albedos which were too high 
to be consistent with previous taxonomic classifications. 

In recent years, Harris (1998) and Delbo & Harris 
(2002) focused on the improvement of the radiometric 
technique for NEAs, resulting in the NEA thermal model 
(NEATM). The NEATM also incorporates a beaming pa- 
rameter ?y, but now it is used as a variable which has to 
be adjusted to produce a fit to the spectral data. Harris 
(JTnUEJ), Harris et al. (fT$9"8"| and Delbo et al. 



3) applied 

the NEATM successfully to thermal N- and Q-band obser- 
vations of a large number of NEAs. Harris & Davies ( 1999 ) 
presented spectrophotometric observations of 3 NEAs and 
discussed the implications for their physical characterisa- 
tion. 

A thermophysical model (TPM) was developed by 
Lagerros 1J1996I 119971 ©98). A number of physical pro- 
cesses were introduced in this model. Instead of the em- 
pirical correction parameters, like the phase angle correc- 
tion in the STM or the beaming parameter in the STM 
and NEATM, the TPM takes the true illumination and 
observing geometry into account and calculates the heat 
conduction into the surface. The non-isotropic emission, 
caused by the distribution of surface slopes, shadows, mul- 
tiple scattering and mutual heating, is described by a 2- 
parameter beaming model with /, the fraction of the sur- 
face covered by craters, and p, the r.m.s. of the surface 
slopes. The TPM, as well as the STM, were mainly used 
for main-belt asteroids. It allows detailed thermophysical 



studies of individual asteroids (e.g., Miiller & Lagerros 
1998)) and predicts for well-known asteroids the thermal 
flux with high accuracy (Miiller & Lagerros f2002fl . Miiller 
( 2002 ) investigated the capabilities and limitations of the 
TPM in the context of NEA thermal observations. 

Here we present N-band photometric and spectro- 
scopic observations (Sect. EJ- The interpretation of our 
2002 NY40 observations was for the first time done with a 
combined NEATM and TPM approach (Sect.EJ). We dis- 
cuss the results of our analysis in Sect. ^ Using the case 
of 2002 NY40 as an example, we also demonstrate the 
possibilities and limitations of applying thermal models 
to asteroids during close Earth ffy-bys. 

2. Observations and Data Reduction 

The 2002 NY40 observations were taken with the TIMMI2 
instrument (Kaufl et al. EH03) at the ESO La Silla 3.6 m 
telescope. 

Table 1. Summary of TIMMI2 observations of asteroid 
2002 NY40 on August 17th/18th, 2002. The observations 
on August 17th were taken between airmass 1.2 and 1.3, 
on August 18th between 1.9 and 2.0. r, A and a are given 
as seen from La Silla. The spectra from August 17th were 
taken with a 1.2" slit, on August 18th with a 3" slit. 



Mid-Time 
No (Day UT) 



Filter 
Band 



[AU] 



A 
[AU] 



1 17 04:48 Ngrism 1.025406 

2 17 04:57 Ngrism 1.025330 

3 17 05:07 Ngrism 1.025245 

4 18 02:02 Nl 1.014650 

5 18 02:08 Ngrism 1.014600 

6 18 02:13 Ngrism 1.014557 

7 18 02:18 Ngrism 1.014515 



0.013973 21.98 

0.013900 22.06 

0.013820 22.16 

0.004535 58.37 

0.004504 58.88 

0.004478 59.30 

0.004452 59.74 



4 http://science.nasa.gov/headlines/y2002/30jul_ny40.htm 
http://www.ing.iac.es/PR/press/ing32002.html 



18 02:23 Ngrism 1.014473 0.004427 60.17 



For all observations we utilized a standard chopping 
and nodding technique to reduce the atmospheric and tele- 
scope background emission. Chop and nod throws were 
10", respectively. For the imaging observations, a pixel 
scale of 0.2" was chosen, and typical on source integration 
times of about 1 minute for each filter were used. Exposure 
times for the low-resolution spectrum (resolution ^160) 
were between 90 and 220 sec. 

The weather conditions were never ideal and parts of 
the two nights were affected by clouds. Table Q] contains 
only the good quality measurements together with the ob- 
serving geometries. 

The photometric Nl-band flux was calibrated against 
measurements of HD 187642 and HD 12929. Colour dif- 
ferences between stars and 2002 NY40 were negligible 
(about 1-3%). For observation No. 4 in Table. ^ we ob- 
tained through standard aperture photometry a calibrated 
8.6 /xm flux density of 22.2±1.6 Jy. 
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TIMMI2 spectroscopic data were reduced and anal- 
ysed using procedures described by Siebenmorgen et al. 
(2004), including an airmass dependent extinction correc- 
tion according to Schiitz &; Sterzik Q2004[l . The absolute 
flux level of TIMMI2 spectroscopic measurements is gen- 
erally better than 5% (Kampgen & Siebenmorgen 2004). 
However, our data are not as reliable mainly due to the 
fast speed of 2002 NY40 in combination with small slit 
sizes. Therefore, we normalised the spectra from August 
18th to the Nl flux density. The spectra from the first day 
were normalised to a model flux (see Sect. El- The spectra 
are shown in Figs, d and [21 The spectral shape is reliable 
in the range between about 8 and 13 /im. 



3. Thermal Modelling 
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Fig. 2. The average of the two best N-band spectra of 
2002 NY40 on August 18th, 2002 (No. 5 and 6 in Tabled), 
both calibrated against HD 187642. The error bar was cal- 
culated in the same way as for Fig.d The residual atmo- 
spheric features were cut out and the spectra were nor- 
malised to the Nl-band photometric flux. 



August 18th took place about 1 hour after the lightcurve 
minimum (J. Licandro, personal communication). 

The constraints from the optical and thermal mea- 
surements can be graphically shown (e.g., Sekiguchi et al. 
2003J with respect to the NEATM assumptions and for- 
mulae (e.g., Delbo & Harris ED02J- Our NEATM analysis 
resulted in an effective diameter of 0.27±0.01km and an 
albedo value of 0.34 + 0.06/ - 0.05 (Fig.©. 



Fig. 1. The best N-band spectrum of 2002 NY40 from 
August 17th, 2002 (No. 3 in Tabled- K was spectroscopi- 
cally calibrated against HD 133774 and against HD 187642 
and then averaged. The 1 er-error bar represents the signal 
processing error together with the 1 cr-average error (on 
the right side). The residual atmospheric features (Ozone 
around 9.58±0.3/im and CO2 around 11.73 and 12.55 /im) 
were cut out. The spectra were normalised to a NEATM 
(j) = 1.0) prediction at 8.6 /xm for the given observing 
epoch. 



We applied the NEATM to determine radiometric di- 
ameter and albedo values from our Nl-band photometric 
flux. For the beaming value we used 77 = 1.5, as recom- 
mended by Delbo et al. (2003) for phase angles larger than 
45°. As asteroid intrinsic input parameters we took an ab- 
solute magnitude of H = 19.6±0.2 mag and a slope param- 
eter G = 0.15 ±0.20, based on an estimated extrapolation 
from 21 lightcurve observations between end of July and 
August 16th, 2002. The H-estimates for August 17th are 
much more accurate, but unfortunately our thermal pho- 
tometry was not reliable for that day. Our observations of 



0.40 

0.35 - 

E 

o 

1 0.30 
n 

b 

LU 

0.25 - 



0.20 



-" "V V "V 


NEATM 0?=1.5): 


\ Optical constraint: 


D, pV determination _ 


\ H=19.6±0.2mag 


- 


\ \. \ 


Thermal constraint: _ 


\ \. 


FD=22.2±1.6Jy _ 


- D,„ = 0.274 ± 0.008 km 


*^<- 


- p v = 0.341 +0.056/-0.049 


'--■■■-■ 



0.1 



0.2 



0.3 
Geom. Albedo pV 



0.4 



0.5 



Fig. 3. NEATM analysis of the Nl-band photometric flux. 
The H-estimate is relatively uncertain and propagates 
mainly into the uncertainty of pv, while the -D e // uncer- 
tainty is dominated by the thermal flux which is relatively 
accurate. The D e ff and pv results are connected to the 
assumption that the model does not introduce a system- 
atic error. 
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Based on these values, we determined a NEATM 
8.6 /im flux density of 5.9 ± 0.8 Jy for the August 17th 
spectra, which we used for normalisation (see Fig.^l. This 
NEATM prediction is based on the Delbo et al. Q2003P rec- 
ommendation to take a value of r\ — 1.0 for data taken at 
phase angles below 45°. 

For both days we calculated spectral energy distribu- 
tions for the mean epochs of the N-band spectra through 
the STM, the Fast Rotating Model (FRM; e.g., Lebofsky 
& Spencer IT^I and the NEATM with -q = 1.0 and 
?7 = 1.5. A nice summary of all formulae can be found in 
Delbo & Harris (2002J). The model predictions are shown 
in Figs. [2 and 

We also applied the TPM by Lagerros ifTU^I H3371 
I1998fl to our observations. We assumed a spherical asteroid 
with an equator-on viewing angle under the given phase 
angle and a prograde rotation with P S id — 19.98 hours 
(Pravec et al. 1277714) . Miiller et al. lfT59"9")> derived / = 0.6 
and p — 0.7 for large regolith-covered main-belt asteroids. 
In addition, wee included a relatively smooth surface with 
a low crater coverage (/ = 0.4, p = 0.4) and a very rough 
surface with 100% crater coverage (/ = 1.0, p — 1.0), 
similar to the studies by Dotto et al. (|2000|l and Miiller 
(2002). For the emissivity, we used a constant, wavelength- 
independent value of 0.9. 

The thermal behaviour of asteroids is strongly coupled 
to the thermal inertia T, in combination with the orien- 
tation of the body in space, its rotation period and the 
sense of rotation. A high T combined with a rapid rotation 
(sub-solar and sub-observer points on the equator) would 
lead to a surface temperature distribution described by 
the FRM. A non-rotating, low T object would correspond 
to the STM temperature distribution. Here, we varied T 
between 10 (regolith covered surface at about 100 K) and 
2500 J rn^s- - 5 !^ 1 (solid rock; JakoskyUJHFJ- 



Table 2. The radiometric diameter and albedo values for 
the different models based on N-band photometry. The 
values are discussed in Sect.0] 



Model 



Diameter [km] Albedo 



Remarks 



STM 


0.20±0.02 


64+ ' 09 


77 = 0.756 


FRM 


0.30±0.02 


29+ 005 

u - zy -0.04 




NEATM 


0.27±0.01 


34+ 006 


77 = 1.5 


NEATM 


0.22±0.01 


S2+ ' 08 


77 = 1.0 


TPM 


0.28±0.03 


0.35±0.02 


r = 1000 


TPM 


0.23±0.03 


0.51±0.03 


r= 100 



In order to obtain a similar diameter and albedo value 
as produced by the NEATM calculations, we had to use 



a high r value of about 1000 J m 



-0.5 



K" 1 . With these 



values we obtained a TPM diameter of 0.28±0.03km and 
an albedo of 0.35±0.02. The errors include the uncertain- 
ties in the measured flux density, the given H-G variations, 
the above specified range of (p,i )-values and a 20 % vari- 
ation in r. 



In Table|2 we summarise the radiometric diameter and 
albedo values for different models applied to observation 
No. 4 (Table HJ). As a combined NEATM/TPM result we 
took D eff = 0.28 ± 0.03 and p v = 0.34 ± 0.06. 
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Fig. 4. N-band spectrum of 2002 NY40 from August 18th, 
2002, with the error bar as in Fig.|5J TPM predictions with 
T = 1000 (upper solid line) and V = 100 Jm~ 2 s^ 05 !^ 1 
(lower solid line) are shown. The dashed and dashed- 
dotted lines indicate the influence of surface roughness 
expressed through p and / for both values of T, respec- 
tively. Default values are / = 0.6 and p = 0.7. 



Figure 0] shows the TPM solution for the August 
18th spectra with T = 1000 J m~ 2 s -0 - 5 K -1 and V = 



lOOJm-'V 



'K" 



For both models the influences of a 



relatively smooth and a very rough surface are shown as 
dashed and dashed-dotted lines. The corresponding figure 
for the August 17th spectrum is very similar. The T = 100 
solution matches nicely, the T = 1000 solution is steeper 
and at long wavelengths significantly above the observed 
spectrum. 



4. Discussion 

Our analysis was guided by the radar size estimate of 
about 250 x420 m (E. Howell, personal communication) 
and the radiometric results (Rivkin et al. 2003) of a 290- 
420 m diameter and an albedo preference of about 0.25. 
All model calculations (see Table 0) produced diameters 
below or at the lower end of these numbers, while the 
albedo was in all cases 0.29 or larger. 

The STM of Lebofsky et al. IPTJ56) has a phase angle 
correction of f3 E = 0.01 mag deg -1 , which is considered to 
be a good approximation for phase angles < 30° (Morrison 
I1977|) . For our photometry at 59° it is therefore not a good 
option. The resulting diameter and albedo values are ex- 
treme and far away from the radar size or the albedo of a 
Q-class asteroid. The STM fit to the spectral slope is as- 
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toundingly good, only for the August 18th data the STM 
prediction falls slightly below the NEATM prediction. 

The FRM led to an acceptable agreement with the 
radar data. But the FRM did not match the spectral shape 
of the two spectroscopic data sets (Figs. ^ and [2J • 

The NEATM (rj = 1.5) and the TPM (r = 1000) pro- 
duced solutions which are comparable to the radar size, 
but again they cannot match the observed SEDs. On the 
other hand, a NEATM with r\ = 1.0 or a TPM with 
r = 100 match the spectra, but produce an unrealisti- 
cally small diameter of about 220-230 m and a very high 
albedo above 0.50. We do not have a convincing explana- 
tion for this, but with only one reliable photometric flux 
it is difficult to find a solution, unless 2002 NY40 turns 
out to be a small, very high albedo object. It might also 
be that strong deviations from a spherical shape seen un- 
der the large phase angle of 58.37° cause thermal effects 
which are difficult to model. On the other hand, one set 
of model parameters fits all spectra, although they were 
taken at very different phase angles. This is true for the 
NEATM and the TPM. 

The NEATM 8.6 /mi prediction for the spectra of the 
first day are 5.9 Jy (r/ = 1.0) and 3.6 Jy (ry = 1.5). The cor- 
responding TPM calculations (taking also the NEATM 
diameter and albedo) give 6.8 Jy (r = 100) and 4.1Jy 
(r = 1000). A reliable photometry for this day would 
therefore indicate which of the models has the better sim- 
ulation of the thermal phase angle effects. 

We cannot confirm a phase angle trend on the beaming 
parameter as proposed by Delbo et al. (2003). A NEATM 
solution with r\ = 1.0 fits the spectra of both days at phase 
angles of about 22° and 59°. 

The T- value is very crucial for our TPM analysis. 



A value of T = 15 J i 



-0.5 



K" 



as it was used for 



regolith-covered main-belt asteroids (Miiller et al. 
leads to a very small diameter and an albedo of about 
0.6. Taking the lunar soil value of V = 39 J m~ 2 s~ - 5 K" 1 
fKeihm ll984[) does not improve the situation much. Only 
large T- values in the order of 1000Jm~ 2 s 



-0.5 



It- 



give 



reasonable results. However, such a large value can only 
be the case if 2002 NY40 has a surface with no or little 
regolith. A dust layer on the surface would always produce 
much lower thermal inertias. 

Through a TPM analysis it is also possible to investi- 
gate the sense of rotation (e.g., Miiller 120021) . We assumed 
a prograde rotator meaning that the terminator is com- 
ing from the cold night side of the asteroid on August 
18th. Taking a retrograde rotator, with the terminator 
coming from the warm day side, would require a T- value 
of about 3000 J m -2 s -a5 K _1 to obtain a diameter similar 
to the radar diameter. Such value would be greater than 
the value for solid rock (Jakosky[1986J and therefore not 
very likely. Again, with only one reliable flux it is difficult 
to obtain a final conclusion, but the TPM analysis points 
to a prograde rotator with a cold terminator after oppo- 
sition. This result will also be useful in case of a possible 
ambiguous solution of its pole orientation in the future. 



The TPM curves in Fig. ^demonstrate that the beam- 
ing model, expressed in p and /, has no big influence on 
the 8 to 13 /mi predictions. Distinguishing between a very 
rough, cratered surface with pronounced thermal beaming 
effect and a relatively smooth surface with no enhanced ra- 
diation at small phase angles was therefore not possible for 
our data set. But the beaming does influence the flux dif- 
ferences between two phase angles (see also Miiller |2*002|) . 
Thermal observations with a good coverage in phase angle 
in combination with a few additional spectra would give 
the possibility to determine T and indicate if the aster- 
oid has a rough, cratered and therefore old surface or a 
relatively smooth, young surface. 

Rivkin et al. (2003) identified olivine and pyroxene in 
the near-IR spectra. The most diagnostic feature in the 
N-band is the Christiansen peak. It is associated with the 
principal molecular vibration band where the refractive 
index changes rapidly. For silicates it appears generally 
as peak between 7.5 and 9.5 /tm (Dotto et al. 2002|. We 
did not see any prominent increase in emittance in this 
range. The feature might have been just outside the re- 
liable wavelength range or the quality of the spectra was 
not sufficient to identify silicate features. 

5. Conclusion 

Guided by the radar size, we found that both the 
NEATM with r) = 1.5 and the TPM with T = 1000 ± 
200 Jm~ 2 s~°- 5 K _1 produce consistent effective diameter 
values of 0.28±0.03km for the time of the N-band photo- 
metric observation. The corresponding geometric albedo 
is 0.34±0.06. 

Neither of the two model solutions give a satisfying 
match to the spectra obtained in two observing sessions, 
separated by slightly more than one rotation period. A 
match of the spectra required either a NEATM with 
rj = 1.0 or a TPM with a thermal inertia in the order 
of 100 J m~ 2 s~ ' 5 K _1 . No phase angle dependence of the 
NEATM 77-parameter can be seen in our data. 

The TPM analysis excluded a dust-covered surface and 
pointed towards a bare rock surface or a very coarse re- 
golith. We obtained reasonable values for the thermal in- 
ertia, the diameter and albedo, but only assuming a pro- 
grade rotation. This means that we saw a cold terminator 
during the observations which took place when 2002 NY40 
was just after opposition, but before its closest Earth- 
approach. 

The example of 2002 NY40 shows the different possi- 
bilities to analyse thermal observations of NEAs through 
the NEATM and the TPM. A carefully planned observ- 
ing campaign under photometric conditions and covering 
a wide range in solar phase angle (without large changes in 
the aspect angles) are ideal for this kind of investigations. 
High quality photometric fluxes are the key to distinguish 
between different model approaches, but a good coverage 
in visual lightcurves is also required to determine H-/G- 
values and the orientation of the object at the time of the 
thermal observations. Additional calibrated N-band spec- 
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tra allow for estimates of the thermal inertia to prove the 
existence of a dust regolith and the sense of rotation. The 
August 17th/18th 2002 fly-by of 2002 NY40 was in princi- 
ple an ideal case for such a study. Only better atmospheric 
conditions for mid-IR observations and a dedicated cam- 
paign in the optical would have provided more accurate 
results. Close fly-bys at distances below 0.01 AU are quite 
frequent (the MPC lists currently 14 such event for the 
period January to May 2004), but most of these events 
are only known on short notice or after the fly- by. In fact, 
there is not even one such event listed the period June 
2004 to December 2009. The next close Earth approach of 
2002 NY40 will be on July 26th, 2005, but at a relatively 
large distance of 0.43 AU. 

Acknowledgements. We would like to thank A. W. Harris 
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NEATM model code. 
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